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The  problem  of  calculating  the  pressure  distribution  on  an  osri'llatiog  plate  in  a 
uniform  parallel  subsonic  she&r  layer  is  idealised  as  a layer  of  fluid  with 
constant  velocity  and  temperature  that  are  different  from  the  free  stream 
conditions.  An  approximate  solution  is  developed  wiui  the  Wiener-Hopf  technique 
for  the  case  that  the  shear  layer  velocity  defect  (displaceaient  thickness/shear 
layer  thickness)  io  small.  The  results  are  valid  .for  arbitrary  length  scale  of 

the  surfape  deflection.  The  problem  is  a model  of  an  airfoil  trailing  edge 
tntsersed  in  a fully  developed  turbulent  layer  and  wake.  The  lowest-order  steady 
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state  solution  is  calculated  for  a linearly  deflected  flap  and  an  exponential 
flap.  Recults  for  she  lintsar  flap  agree  with  the  known  solution  for  i finite 

length  air-foil.  The  lowe«t-order  unsteady  exponential  flap  and  the  first-order 

steady  state  exponential  flap  are  also  solved  in  detail.  The  first-order 
pressure  decreases  as  the  s^juare  root  of  the  distance  fro®  the  trailing  edge 
like  the  lowest-oreder  result,  and  in  accordance,  with  a full  application  of  the 
Kutta  condition.  The  aairlicude  of  the  correction  is  of  order  t n 1/  where 

£.  the  ratio  of  the  shear  layer  thickness  to  flap  length.  Xt  is  shown  that 

the  total  lift  on  any  section  of  che  trailing  edge  of  length  L and  flap  length 
i increases  like  (W^  as  L increases  v<tho»:t  bound  in  the  absence  of  a 
sb^r  layer.  An  experinent  to  detenu ne  how  the  shear  layer  lisilts  the  growth 
of  lift  is  suggested. 
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I.  INTRODUCTION 


In  a previous  long-range  effort  (see  Ref.  1,2)  under  Air 
Roree  sponsorship,  the  author  has  developed  aii  integral  approach 
for  Including  the  effect  of  the  boundary  layer  in  the  convent Icnal 
unsteady  aerodynajnlc  analysis.  The  original  objective  of  the  pro- 
=gi‘<aia;  to  develop  the  theory  to  a point  where  it  could  be  used 
in-  the  analysis  of  supersonic  panel  flutter.  This  objective  has 
been  accomplished.  The  Integral  technique  has  been  sv.ecessfuliy 
Integrated  into  a pariel  flutter  computer  program  (Refs.  3,^,5) 
snder  Air  Force  and  NASA  sponsorship.  The  effect  of  the  boundary 
layer  on  panel  flutter  has  been  calculated  and  compared  with 
etperls»ntal  results.  The  agreement  Is  remarJcably  good. 

Because  of  the  success  of  the  Integral  approach  in  solving 
the  panel  flutter  problem,  a shorter  range  program  ms  Initiated 
to  investigate  related  applications  of  the  theory.  The  particular 
probls.%  of  noise  generation  and/or  reflection  by  a panel  surface 
in  a shear  layer  was  a possible  source  of  application.  The  trailing 
edge  problem  was  a second  possibility.  Our  work  on  the  acoustic 
problem  has  resulted  in  a revolutionary  new  theory  of  aerodynamic 
s?^nd  generation.  This  work  is  reported  in  two  published  documents 
(see  Sefs.  0,7),  Our  final  report  is  concerned  solely  with  the 
trailing  edge  problem. 

To  the  author’s  knowledge,  the  first  work  on  the  appllcatJof: 
of  shear  layer  aerodynamic  theory  to  lifting  surface  problems  is 
due  to  Sowell  and  Ventres  (Ref,  8).  In  a subsequent  study.  Ventres 
{Ref.  f)  developed  the  shear  layer  kernel  function  in  detsil  for 
steady  incompressible  flow  and  calculated  the  effect  o-f  the  shear 
layer  on  the  lift  curve  slope  and  center  of  pressure  of  various  two- 
and  three-dimensional  airfoils.  Ke  showed  that  the  .lift  curve  slcpe 
is  decreased  while  the  center  of  pressure  is  unaffected  by  the  shear 
layer.  In  a related  study,  Williams  (Ref.  JC)  .has  developed  a 
general  approach  for  solving  a generalised  class  of  singular  integrj’ 
equations  that  result  when  the  shear  layer  is  introduced  into  the 
usual  aerodynamic  theory.  More  recently,  Dowell,  Chi  and  Willla.T.t 


il)  have  been  atteopting  to  extend  the  kernel  function  con- 
cept to  unsteady  flow  by  expanding  for  sisall  frequency  about  the 
steady-state  solution  of  Ventres. 

The  present  work  is  closely  related  to  the  work  of  Dowell 
and  his  co-workers  in  that  we  start  with  the  sane  basic  TOdel  of 
the  shear  layer  and  focus  on  the  solution  of  a lifting  surface 
.crdblem  as  opposed  to  the  panel  aerodynaaiic  problem.  Our  approach 
is  soisewhat  different,  however,  in  that  we  attempt  the  development 
of  an  analytic  solution  for  a very  special  airfoil  and  shear  layer 
modei.  Specifically,  we  treat  a semi-infinite  airfoil  and  replace 
the  actual  shear  layer  by  a flow  with  constant  'velocity  and  teiq>era- 
ture  that  are  different  from  the  free  stream  conditions.  Ah  approxi- 
mate solution  is^  developed  with  the  Viener-Hopf  technique  in  the 
limit  of  small  shear  layer  velocity  defect.  Explicit  results  are 
given  for  the  case  of  a linearly  deflected  flap  and  an  exponential 
flap.  The  analytic  results  should  be  xiseful  for  evaluating  our 
Integral  approach  and  other  approximate  schemes  for  solving  kernel 
function  models. 
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II..  (FORMULATION  'OF  THE  'TRAILING  FDGE  FE03LEM 


.Siya?temeiit  ■&!  ^'he  Fr-obleni 

Me  'OG/nsider  tlhe  pr.oblem  iii;ust>ra1>ed  in  Fig..  1..  .A  -two- 
-dimensional ilat  iplafce  “'xt.ends  .along  tlie  ;ne,gati-ve  atxis  witlh  a 
inaiTing  -edge  ;at  fche  .origin..  '.The  .plate  .undergoes  .small  oscillations, 
Tihe  flo.w  .is  lunii.orm  ^subsonic  lexcept  in  a tLin  layer  .oi  tlhiciknes-s  6 
near  ihe  :x-.axis  where  the  imean  w.elocity  is  -a  function  of  y only,. 
The  ibasic  iproblem  is  to  (determine  the  .pressure  distribution  on  the 


'The  model  -jprohlem  w.e  hawe  adopted  .i;s  an  idealization  <.of  the 
trailing  adge  nf  .an  -.os.ci Hating  aii’foii..  The  primary  ;aim  is  to 
determine  the  affect  -of  the  ahear  .layer  nn  the  surface  pressure 
near  the  trailing  adge,.  The  .solution  of  the  corresponding  inviscid 
problem  is  iknown  ((‘Ref..  121)  and  .is  readily  obtained  with  the  iWiener- 
•Hopf  technique,.  In  the  present  work  we  .show  how  the  same  technique 
■can  he  .used  to  obtain  an  approximate  .solution  for  a simple  two-layer 
‘.model  of  the  shear  layer.. 

Basic  Equations 

The  -.equations  -for  the  perturbation  pressure  p and  -vertical 
■?/elocity  -.w  are  v(,Ref,.  1)  in  dimensionai  form:: 


^ - di-v  (0  grad  p)  = 2^'  4^ 

,-af  .Bt^ 


v;here 


Sr  + '0  = 0 

Bt  -3y 


pert urb a t ion  pre  s s ure/dens i ty 
free  stream  sound  speed 


■V  ^ + ,5— 
:3x  3t 


( 2 . ? ) 


3 


ir  = 17 (y3  .shear  layer  -velocity  profile 

T„lyD 

nonnall.sed  temperature  profile  in  the 
shear  layer 


0 * 


- o 


(2.3i 


l!?e  assiune  the  motion  of  the  plate  to  be  simple  harmonic  in 
time  (q  -*■  so  that  p and  -m  are  singple  harmonic  at  every 

point  in  space,  Purthermore,  -we  Pourier  analyze  the  x dependence 
of  eaeh  dependent  variable.  Thus, 


and 


(2ir) 


T72 


q<xO 


(2ir0 


IT? 


f 


dx  qixO 


ix  e q(a) 


12M 


■where  aj  is  either  p or  w , The  equations  for  the  Pourier 
amplitudes  become 

iC.a¥  w)  * er^p 


^ ^ , (at  + J£l  - 
dy  0 


■where 


^2  ^ (at  uV 


afe 


(2,53 


(2.bl 


Boundary  Conditions 

In  the  region  outside  the  shear  layer  (see  Pig.  1)  ve  obtain 
a single  equation  for  p ; i,e,. 


H-rip^o 

dy 


(2,7) 


■wnere 


r 2 


r * 


7 


(2,fi) 


I 

I 

I 

I 

I 


¥e  require  tnat  the  pressure  either  decays  exponentially  or 
corresponds  to  outgoing  (acoustic)  va^?es^  Thus, 


p - for  y -*  a 


(^.5) 


and  the  boundary  condition  is  satisfied  ±f  we  choose  the  branch  cuts 
for  the  conplex  function  as  shown  In  2.  The  branch 

points  a,  and  =Lr-e 


to 


= + 11 

1 

OC  Q3 


to 


^ °2  "a  + 

St  ao 


and 


r„  = S„(o  - 


= U . 


(2,10! 

(2.11) 

(2.12) 


For  reasons  that  will  become  clear  later,  we  assume  that  lo  has  a 
small  positive  complex  part  so  that  is  sli^tly  above  the  real 

axis  and  is  slightly  below  the  real  axis  as  indicated  in  Fig.  2 

On  the  real  axis  we  can  write 

r = S a-a.  ja  — a,,  a>a, 

OD  00  ^ i.  * * C * 1 

oo 

or  a < 


= “ is  la  - a.. 


,1/2,  00,7^2 

1 - Qul 


a,  < a < 


(2.15; 


which  results  satisfy  the  required  boundary  condition  for  y -*  ® . 

Additional  boundary  conditions  in  the  shear  layer  are  the 
following.  The  pressure  must  be  continuous  at  all  points.  If  the 
mean  flow  properties  are  continuous,  then  the  velocity  w is  also 
continuous.  If  we  admit  a discontinuity  in  v , then  continuity 
of  particle  position  requires  that 

(aV^  + to)w  = (rxV  + to)w^  (2.1^) 

where  the  + and  - subscripts  refer  to  the  points  slightly  afcove 
or  below  the  y location  of  the  discontinuity. 


6 


&e  bouncary  conditions  on  the  x axis  are  as  follows: 


w{x,0)  » ¥(x)  * II  + iwf  for  X < 0 (2.15) 

p(x,0)  * 0 for  X > 0 (2.16) 


where  f(x)  is  the  deflection  of  the  plate.  Finally,  we  require 
that  the  pressure  tend  to  zero  at  the  trailing  edge  of  the  plate. 
Thus,  we  have  invoked  a fonn  of  che  sCutta  condition. 

Coffiiaents 

solution  of  the  problen  thus  f*orjsulated  is  difficult  for 
an  arbitrao’y  shear  layer  profile,  .Mach  number,  frequency,  etc.  Thus, 
we  consider  a sln?>lifled  probleai  by  choosing  a two-step  shear  profile 
as  rhpwn  in  the  next  section.  ¥e  remark  that  other  slnpliflcations 
om/also  be  made.  For  exan?>le,  in  the  Incompressible  case  we  have 

0*1 

- a~/af  <2.17) 

and  from  (2.5)  we  find  that 


ds 


* 0 


(2.18) 


For  a linear  velocity  profile,  the  equation  for  w can  be  solved 
exactly  in  terms  of  exponential  function  and  the  pressure  follows 
by  Integration.  Recently,  Goldstein  (Ref.  13)  has  considered  the 
compressible  problem  Tor  which  there  is  also  an  exact  solution  for 
a linear  velocity  profile. 
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III.  SOLUTION  0?  THE  TWO-LAYER  MODEL 


The  Reduced  Problem 

We  consider  the  simplified  shear  layer  shown  in  Fig.  3.  The 
"Velocity  and  temperature  are  assumed  to  be  constant  in  a layer  of 
thickness  6 . Free  stream  conditions  hold  for  y > 6 . 


where 


In  the  layer  the  pressure  p satisfies  the  ordinary 
differential  equation. 

lx  - r^p  » 0 
dy^" 

r 

^ _ 1.2  (aV  + ur 

‘ L 

2 2 

a'^  * a^e  speed  of  sound  in  the  layer 
The  solution  of  (3*1)  is 

p * Be"^*^  + Ce^^'  0 < y < £ 


(3.1) 

(3.2) 

(3.3) 

(3.^) 


where  the  definition  of  F in  the  complex  a-plane  is  precisely 

the  sa.me  as  f'’.r  F (see  Pig.  2).  We  have 

00 


r 

B 


6(o  - 


(1  - m2)1/2 


1/2 


M - V/a 


a 


i 


*^2 


(3.5) 


The  velocity  in  the  layer  follows  by  diff entiation  of  (3.^)  and 
use  of  (2.5^;  i.e.. 


9 


iSi 


i?  = ierCBe"^^  - Ce^y)/(aV  + w) 
For  y > 5 , we  have  (see  (2.9)) 


) < y < 6 


(3.6) 


-r^y 


p * Ae 


iA  e 


w * 


V + 


The  pressure  and  pai*ticle  position  must  continuous  at  y 
(riote  (2.1^).)  We  thus  solve  for  B and  C in  terms  of  A 
finally  calculate  the  ratio 


where 


C - r - „«-2r6 
g = E - oe 


rOCaV.+u)  - r^(aV+u) 

J- 

r0(aV„+u)‘^  + r (aV+u) 


For  y » 0 , we  get  from  (3.^)  and  (3.6) 

p (1  + I)B(a) 

. i (txV  t_a>2  w . r(l  - 1)8(0) 

we  formally  invert  this  pair  of  equations  to  obtain 


p(x) 


1 ^ 
^ Dx 


(2T7) 


IT? 


(211) 


172 


«D 

/- 


(1+2)  • B{a)  • do 


(1-2)  • 3(o)  • do 


(3.7) 

6 . 
and 


' I?  • ^ * 


(3.9) 


(3.10 


(3.1) 


10. 


wnere 


D 

m 


V §3^  + i« 


(3.i: 
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These  are  the  basic  equations  that  we  use  to  study  the 
trailing  edge  problem.  The  function  i3<a)  is  the  essential 
unknown.  We  have  cast  the  problem  in  the  above  form  to  facilitate 
application  ef  the  Wiener-Hopf  technique.  In  the  complete  absence 
of  a shear  layer  we  have  I » 0 and  Wiener-Hopf  solves  the  problem 
exa^ctly  as  we  shall  see  below. 

Application  of  Wiener-Hopf 

We  apply  the  boundary  conditions  on  the  x-axls  (see  (2.15)  and 
(2>l6))  to  obtain  the  following  pair  of  dual  integral  equations: 


g-io(X(i  ^ £)B(a)da 


0 , 


X > 0 


(3.1^) 


-lax 


(1  - I)B(a)da 


1 DW 

S' 57 


X < 0 


(3.15) 


Next,  we  follow  the  Jones  method  of  applying  the  Wiener-Hopf 
technique  (see  Ref,  12)  to  partially  symmetrize  the  integrands  in 
(3.1^)  and  (3.15).  First,  we  replace  x by  x 5 in  (3.1^)» 
multiply  by  ”scme  function"  N^(^)  and  integrate  from  0 to  ® . 
We  get 


1 

(2x) 


175” 


g-iaxjj  ^ £)B(a)da  » 0 


X > 0 (3.15) 


where 


N (a) 


e-laxjj^(x)dx 


o 

and  N_(a)  must  be  analytic  in  "some"  lower  half  of  the  complex 
plane.  Slmilanrly,  we  replace  x by  x - 5 in  (3.15),  multiply 
by  "some  function"  N2(C)  and  integrate  to  obtain 


t “i 


12 


1 


(2w)- 


)r(l  - r)B(a)da 


X < 0 


where 


00 

i^(a)  ‘J  e^®*N^ 


(x)dx 


and  N^(a)  is  analytic  In  ’’some"  upper  half  of  the  complex 
plane . 

Next,  we  construct  the  functions  N^(x)  and  N^Cx)  by 


choosing 


N (a)  = l/(a  - o^)‘ 


N^(o)  « 1/(0  - (3.2C 

where  o^^  and  Og  are  the  branch  points  of  the  function  r(o) 

(see  (3.5)  and  Pig.  2).  Note  that  N_(o)  Is  analytic  in  the 
lower  half  plane  (Im  a < Im  o^)  and  N^(o)  Is  analytic  in  the 
upper  half  plane  (Im  a > Im  Og)*  It  is  essential  for  the  analytic 
continuation  used  in  the  Wiener-Hopf  technique  that  these  two  half 
planes  overlap.  This  was  the  reason  for  assuming  a small  positive 
complex  part  of  u . We  also  remark  that  the  particular  choice  for 
N_(a)  and  N_j^(a)  establishes  the  cnaracter  of  the  pressure  near 
the  trailing  edge.  We  shall  see  that  p(x)  tends  to  zero  as 
for  X -►  0“  in  accordance  with  the  Kutta  condition. 

With  (3.20)  the  functions  N^(x)  and  N2(x)  are  easily 
evaluated  (see  Appendix  A).  We  get 

io,x+l7t/ii 

Nj^(x)  ■ (irx)"  ' e , x > 0 

. yp  “iOpX**i  it/4 

Np(x)  » (xx)"-^''^  e , X > 0 (3.21 
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Both  N,  and  N-  are  zero  for  x < 0 . 

^ la,x 

Finally,  we  multiply  (3*16)  by  e ^ and  differentiate 

with  respect  to  x . Our  dual  Intesi'al  equations  (3»1^)  and  (3*i5) 

are  reduced  to  the  following: 


(2it) 


1/2 


(2it) 


m 


00 

/ 

*00 

00 

/ 


g-lax^j^  + r)Q(a)<3a  = 0 , x > 0 


e-iax(i  ^ £)Q(a)da 


where 


00 

o 


Q(a)  * (a  - aT)^^^B(a) 


X < 0 


Once  Q(a)  is  known,  the  wall  pressure  is  calculated  with 

1 


p(x) 


(2t») 


T72 


/ 


e-lo*  d« 

(a  - 


Exact  Inviscid  Solution 


(3-22) 


(3.23) 


(3. 25' 


In  the  absence  of  any  shear  defect,  we  have  1=0  and  the 
pair  (3*22)  and  (3.23)  yield  the  exact  solution  of  the  problem;  i.e. 

o 


Q-(a) 


o 


(2x) 


T7? 


/ 


^iax  , 1 

e ax 


/ 


o 


N,(C) 


dC 


f ^ 


where  B * 8 and  0=1.  For  the  surface  pressure  we  ge* 


oo 

/ 


-lax 


Qo«.! 


(a  - 


da  , X < 0 (3. 


1^1 


We  can  reduce  these  results  further  with  the  convolutiort 
integral.  For  (3.27)  we  get 


Po(r)  » I - x)Q^(4)d5  , 


X < 0 


■fhere 


° (2isy^ 


e“^®^Q^(a)da  , C < 0 


Now  substitute  (3.26)  into  (3.29)  and  invert  the  or^er  inte- 
gration to  get 


cc 


X ' u 


Thus,  the  exact  pressure  on  the  surface  becofnes 


J^N^(s  - x)d3  J* 


where  and  N2  are  defined  by  (3.21). 

For  any  specific  downwash  distribution,  W(x'  , the  ralrulation 

of  p^(x)  has  been  reduced  to  a double  ouadx at --r.; . *r 
o 

specific  examples  in  the  subsequent  section. 

An  Integral  Equation  for  Q(a) 

The  obstacle  that  prohibits  the  exact  solution  ;or  iv'x  in 
•;3.23)  and  (3.23)  is  the  complex  function  Ka)  . We  ' ar  dlspujse 
our  difficulty  somevxhat  by  transposing  the  two  terms  thir  : r.»  ain  l 
to  the  right-hand  side  of  the  equation  and  rorna:ly  irve:  ♦ t nf'.  The 
result  is  a singular  Integral  equation  for  the  jnkrrwi.  .. .a»  ; l.a*  , 


QCa)  = Q^(a) 


■t/ 


Z(o ' )Ql a ^ * 


I 

i 

I 


where  Qq(«)  is  given  by  (3*30)  and  the  integral  is  a standard 
Hilbert  transform  (Cauchy  principal  value  integral). 

Since  we  still  do  not  know  Q(a)  , we  have  only  reduced  two 
integral  equations  to  one  with  the  last  trick.  However,  the  form 
of  (3*28)  suggests  that  we  might  iterate  for  Q . If  £(a)  is  in 
some  sense  small,  we  can  expect  to  obtain  useful  results.  We  post- 
pone the  Investigation  of  ^ for  a moment  and  formally  derive  a 
first-order  result.  We  have  to  first  order  in  £ 


- Qo  + Q' 


and 


1 C 

' a'  - . 


do' 


*5 


II 


I 


\ 


Next,  we  evaluate  the  first-order  pressure  with  (3.25).  We  have 


(1  + Z)Q^  • (1  ♦ £)Q^  -J-  Q'  ^ OUn 


(•5 


and  we  write  the  nressure  in  the  form  of  a convolution  integral, 
Thus, 


wnere 


Consider 


P(C)N^(e  - x)dC 


P(x) 


X 

00 

J + £)«„  + 


(3.35? 


(3.36: 


?'(x) 


sgn( 

(2ff) 


OC 

I 

ao 

I 

00 

n(-x)  I 

* 


e”^®^Q'(a)da 


I(a')Q^(a')da' 


00 

i?r  J a - a 


do 


^“^£(o)Q(a)da 


( 3.3?' 
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It  follows  by  inspection  of  (3-^3)  that 


-1  , o(a^  _)  = +1 


.a(-u)/V  .)  = -i  ^ o(-u)/V)  * +1 

ac 

Por  a -►  ± •“  the  paraTnete.r  o becomes  independent  of  a 


(3.50) 


2 2 
60i^„  - e V"" 

flD  flO 


B0V  + B V 
00  00 


-2S6jal 


for  lal  - » 


(3.51) 


I ~ oe~-““'“'  (3.52) 

Ve  also  remark  that  for  the  steady-state  case  the  expression 
(3.51.)  is  valid  for  all  a . 

The  function  I is  real  for  a > a*  or  a < . It  is 

bounded  uniformly  on  the  real  axis  and  never  exceeds  unity  in 
absolute  value  (see  Pig.  ^).  In  the  steady  state  it  is  hounded  hy 
o . Por  low  Mach  number,  we  have  approximately 


1 - v^/v^ 

OD 


U.  + V^/V" 


1 - 3L  = -A  _ 

^ V N + 1 


(3.53) 


(3.54) 


vhere  6 is  the  displacement  thickness  and  N is  the  index  of  a 
power  law  \'elocity  profile.  Thus 


)L  - 1 1 W 

V "-^"n  + i'n  + i 
00 


2N  + 1 


2N‘  + 2N  + 1 


.13  for  N » 7 


(3.55) 


(3.56 
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Sketch  or  Complex  a Function 


For  a steady-state  deflection  a is  of  the  order  of  the  boundary 
layer  velocity  defect,  and  we  expect  the  iteration  scheme  to  work. 

A rigorous  proof  could  be  developed. 

Since  I is  only  bounded  by  unity  in  the  unsteady  case,  we 
can  not  prove  a priori  that  the  iteration  scheme  will  work.  However, 
it  seems  plausible  that  it  will,  at  least  for  sufficiently  low 
frequency.  The  question  is  whether  we  can  ultimately  calculate 
results  for  large  enough  frequency  to  see  any  significant  nonsteady 
behavior  near  the  trailing  edge.  We  remark  finally  that  when  1=1, 
the  two-step  shear  layer  has  mstable  elgensolutions  and  when  Z = - 1 
the  wake  has  unstable  elgensolutions.  This  fact  may  present  some 
coiipli cations  in  obtaining  a unique  solution  of  the  unsteady  lifting 
surface  problem.  Further  analysis  of  this  important  point  should  be 
carried  out  before  we  Imbark  on  a numerical  solution  of  the  unsteady 
problem. 


IV.  APPLICATIONS 


We  Illustrate  the  foregoing  theory  with  two  examples.  First, 
we  calculate  the  lowest-order  solution  for  the  linear  flap  at  zero 
frequency.  Second,  we  consi<3er  an  oscillating  exponential  flap  for 
which  we  calculate  lowest-  and  first-order  results. 


The  Steady  State  Linear  Flap 

Consider  the  steady  state  deflected  flap  shown  in  Fig.  5.  The 
surface  deflection  is  given  ty 

f(x)  = - 4)(x  - Xj^)H(x  - x^)  (^.1) 

where  H(x)  is  the  Heaviside  step  function.  The  downwash  and 
surface  acceleration  are 


W * Vf'  * -♦VHCx  - x^) 


(^.2) 


and 


. VW-  . -*V=5(x  - x^) 


(4.3) 

Where  6(x)  is  the  Dirac,  delta  function.  It  is  understood  that  x 
is  negative  in  all  of  the  subsequent  formulas. 

The  lowest-order  surface  pressure  is  obtained  with  (3.31) • We 


nave 


o ^ 

*X  *0 


o * 

- |-Q-  J*  N^(s  - x)ds  J*  -5-Xj^)d4 


it,T 

AS 


o 

; 


K^(s  - x)N2(s  - Xjj)ds 


Max  ( X , 

.o 


£V“ 

II  se 


/ 


/s  - X /3  - Xi 


Max(x.xj^) 

■ 'S®  " L lx  - XV  l^-  - J 


With  B = S>.  and  6 


the  last  result  is  the  exact  solution  of 


the  linear  deflected  flap  in  a nonshear  flow.  The  pressure  is 
plotted  in  Fig,  6. 


Remarks 

We  note  several  asynptotic  results.  First,  the  pressure  tends 
1/2 

to  zero  as  (x/Xj^)  ' at  the  trailing  edge  and  has  a logrithmic 

singularity  at  the  hinge  line  in  accordance  with  known  results  for 

the  finite  airfoil  (see  Ref.  13*  e.g.).  Also,  the  pressure  tends 

1/2 

to  zero  as  (xjj/x)  for  x . This  means  that  the  total  lift 

on  the  plate  is  noncalculable.  In  fact,  the  lift  on  any  length  L 
of  the  trailing  edge  tends  to  infinity  as  for  L -*■  « . This 

is  in  agreement  with  the  known  exact  solution  (Ref.  1^,'  for  the 
lift  on  a finite  airfoil  with  chord  c at  zero  angle  of  attack,  but 
with  a finite  flap  deflection.  The  total  lift  tends  to  infinity  as 
/c  . We  shall  discuss  the  implication  of  this  singular  result  in 
more  detail  after  we  consider  the  next  example. 

The  Exponential  Flap  (Zero  Order) 

The  linear  flap  is  a simple  example  of  the  lowest-order  steady 
state  solution.  However,  the  algebra  becomes  very  tedious  when 
higher-order  or  even  unsteady  results  are  sought.  For  this  reason, 
we  consider  the  unsteady  exponential  flap  illustrated  in  Fig.  7.  The 
flap  deflection  mode  is  assumed  to  be  exponential;  i.e.. 


sc  that 


f(x)  » - 


W « ^ = -aV(l  - ik)e*'^ 


V-i.5) 


wnere 


§1  « - I (1  - ik)^e*'^^ 


k » u£/V  , reduced  frequency 


and  ^ is  an  effective  flap  angle;  i.e.  $ * -f'(0)  . 


Distribution  on  th«  Linear  Flap 


,^_  ' - - - - ' --^  ■ ■ 


Our  next  step  is  to  calculate  the  function  Qq(J^)  with  (, 


and  Q^(o)  with  (3.26).  We  have 


w 

o 

- *V^(l-i^)^  I e' 

■ " 30i  J 


e/t 


With  N^(x)  given  by  (3.21),  we  get 


00  00 

j e“^^^N2U)dC  “ J*  e 


-i/i  e 


-ittpC-i^/^  1/2  -ix/4 


dC 


1 ' e 


(1  + ia2l) 


IT? 


so  that 


Q*(x)  * (2x)^^^  I 

o *> 


where 


The 


B0(l+ia2fc)^'^^(2x)^^‘^ 
comolex  function  Q„(tt)  follows  from  (3.26): 


o 

I 


«o'“’  ” 


■if 


lax  x/i 


e ' 'dx 


C 

1 + ioT 


Thus,  is  a simple  pole  at  a » i/s.  . 

The  lowes-  'der  pressure  follows  from  (3.28).  We  get 


3.29) 

(i4.8) 

(4.9) 

(4.10) 

(4.11) 


(4.12) 
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P„(x) 


U - x)Q^(t)<Jt 


-x/i 

. p /7  U/^  x/l  r 
= C/^e  ‘eje  -~ 

o 


The  integral  in  (4.13)  can  be  expressed  in  several  different  forms 
with  the  complex  error  function.  However,  the  simplest  form  for 
computational  purposes  is  to  write  the  real  and  imacinery  parts  out 
separately.  Thus,  we  get  after  substituting  for  C : 


where 


G(z)  * e 


ReG(z)  * e 


•'o 

■’7’  •' 


l+iaj^£)t‘ 


ImG( 


3,  - 


e cos  (c»,£t  )dt 


2 

e sin  (ci^)lt^)dt 


(4.14 


(4.15 


(4.16 


The  oscillating  exponential  flap  can  thus  be  solved  comi.,'etely  in 
lowest  order  and  expressed  in  terms  of  simple  functions. 

It  is  interesting  to  consider  the  steady  state  exponential 
flap . We  have 


Po<  = > ' - # p[c-»/«-'^] 


Where 


P(z)  * e 


2 r 2 


(4.17 


(4.18 


is  Dawson’s  integral  (see  Ref.  15).  The  pressure  profile  corres- 
ponding to  (4.18)  is  plotted  in  Pig,  8.  The  pressure  tends  to  zero 


as  as  X ♦ O"  9 and  tends  to  zero  as  /-i/x  for 

X . Thus,  the  total  lift  is  again  infinite,  a result  that 

does  not  seem  to  depend  upon  the  details  of  the  surface  deflection. 
The  pressure  has  a sninimum  as  indicated  In  Fig.  8.  Thus 

2 

l«ln  Pol  * -SW  /||f- 


I = - .85^11 


(^.IS) 


First-Order  Exponential  Flap  (Steady  State) 

Because  of  the  simplicity  of  the  function  Qq(o)  > we  can 
derive  first-order  results  for  the  exponential  flap  in  steady  state. 
In  principal,  the  unsteady  first-order  result  can  be  calculated  but 
the  details  are  extremely  tedious.  The  pertinent  formulas  of 
Section  III  are  summarized  belowi 


Pi(x)  » Pq(x)  + p'(x) 


p'(x) 


P'(x) 


r*(x) 


J*  P'(UN^ 

— o> 


- x)de 


Q*(Oz\x  - C)dC 


(^.20) 


wnere 


Z(a)  = 


V 01  = 0 


(4.2j 


2 2 
sev^  - sjr 

O)  CO 


B0V„  + 8 

<jo 


Again,  it  is  understood  that  x < d in  all  of  the  formulas 
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The  -reason  itor  Tthe  ;simpli-cit;y  ‘Of  nthe  ;steacly  .stat-e  <case  is 
Tthat  <0  ^ defined  ib;y  is  independent  <of  >a  'Thus^  -.we  -can 

^evaliuate  each  nf  ihe  functions  needed  in  the  .series  of  convolution 
int'pgrais  in  ((■^..;2,0h.  iFirst^  me  caic-uiate 


t*lx)  ^ ^ 


“iax  ;2B6i|tajl 


-,23'6ct 

•e  <cos  'OX 


((2^5 


■IJo.B'fi 


IWe  'Substitute  this  result  together  with  into  the  third 

formula  '.of  (('i|...20:)  to  Obtain 


■8C.3a'6 


e^^^  d£ 


+ ilB  -6' 


8CoB6 


e~®ds 


:S  + i<:B  6 /£' 


The  last  result  can  ,bve  expressed  in  terms  of  exponential  integrals, 
but  it  is  more  convenient  to  leave  it  in  the  simple  integral  form 
for  the  moment. 

Pinally,  we  substitute  (^.2^)  together  with  the  definition  of 
(see  (3.2?.))  into  the  expression  for  p'(x)  . We  get 


P'(x) 


P'(^)N, (e  - x)dC 


8CoB6 


J (C  - x)' 

X 


e”''  ds 


j.  /£" 


(^.25) 
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or  after  eome  manipulaciori  cf  the  double  integral 

,2  r.  ..a/2 


U,  =-S-/|  3.p[(-f)^  ,ifi] 


v;here 


y(ZjE)  = e 


-z 


7 « • r / 


t^-Z^ 


ds 

T 

e 


The  final  expression  for  the  first-order  pressure  is 


(4.26) 


(i|..27) 


(i|.28) 


The  function  ?(Zje)  is  a generalization  of  Dawson’s  integral;  i.e.. 


i;’ 


Lim  F(z,e)  * F(z) 

E-“0 


(^.29) 


V7here  ?(z)  is  given  by  (4a8).  Thus  we  can  also  write  (^..28)  in 
the  form 

r2 


where 


p (x)  = - -^^^7(1  * 2o)P 


F*(z 


l-rj 

-2oP*[(-  i)"^^  2|ijj 

.<%.  i / ^ 


-s 


ds 


(^(.30) 


( . 31 ) 


t^-z^ 


mm 


fl 


tf 

mm 


Asymptotic  Results 

There  are  two  independent  parameters  in  our  first-order  shear 

layer  r-esu1t.  The  basic  expansion  parameter  is  a and  terms  of 
2 

order  o nave  been  neglected.  This  parameter  is  truly  a property 
of  the  shear  layer  as  we  have  shown  in  Section  III  (see  (3-53)  to 
(3«56)).  The  second  parameter  is  2B6/S.  and  it  is  the  ratio  of  the 
shear  layer  thickness  to  the  length  scale  of  the  flap.  We  can  vary 


iC-% 
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t-bis  parameter  independently  of  o . For  example,  when  6/SL  ■*  0 
the  pressure  is  given  by  the  first  term  in  (^.30);  i.e-. 


Lim  p^(x)  = - 
6/H-O 


(1  + 2o)F 


[(-  rf"] 


In  the  complete  absence  of  a shear  layer,  we  have 


2oY^ 

p (x)  = ^ F 

° /i-6 


[(-  i)‘"] 


SO  t-hs.'fc 


y = ;r-rrr  ” 2o) 


P„^x, 


(^.3^ 


The  importance  of  the  ; actor  2o  is  now  clear.  Without  this  term 
(ii.32)  would  tell  us  that  the  lowest-order  pressure  is  to  be 
obtained  by  replacing  the  entire  flow  by  a uniform  flow  with  the 
velocity  and  temperature  (or  density)  of  '^’he  shear  layer.  This  is 
a gross  over-correction  of  the  inviscid  result  and  roust  not  be  used. 
The  correct  factor  to  use  is  y given  by  (^.3^). 

2 

For  low  Mach  number  M <<  1 we  can  exoress  y in  terms  of 

QD  * 

the  shear  layer  velocity  defect;  i^e.. 


(^.35; 


Then 


1 + 2 


, . 4 r. . ,ri^)] 

< [ \v;  + vv  j 

= (1  - i)^fi  + ■■  r ^ I 
L 1 - + A /2J 


= [1  - 2A  + 0(A“)][1  + 2A  + GU")] 
= 1 - O(A^) 


The  correction  due  to  the  shear  layer  is  0(A*',  and  is  therefore 
r._-r.talculable  with  first-order  theory  when  {.  L - 0 . 


For  6/£  * 0(1)  we  must  in  general  use  (-4.30)  to  calculate 
the  shear  layer  correction.  For*  small  6/ 1 we  can  estimate  the 
order  of  the  correction  with  the  asymptotic  results  derived  in 
Appendix  B.  We  have 


-if  9^. 

ao 


e^^GCn.e) 


(4.37) 


1 - - tan~^  - 
IT  n 


- - ie)  - e^^E^(-Ti  + ie^ 


(*4,36) 


where 


|arg  z]  < ir 


(4,39) 


For  e -►  0 we  have  further 


3-0  J I L 


(h^+e^)- 


E^(ti) 


where 


E^(n) 


OD 

-f-? 


(4.40) 


(4.41) 


In  the  vicinity  of  the  trailing  edge,  we  get 


F*(^z,e)  * ^ ^In  i + 1 - ^ /z 


(4.42) 


E -►  0 


where  y * ,57721  is  Euler’s  constant.  Thus,  the  first-order 
pressure  tends  to  zero  as  near  the  trailing  edge.  However, 
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the  magnitude  of  the  correction  is  0(e  In  1/e)  where 
e = 2S6/S.  . We  conclude  that  any  expansion  of  the  shear  layer 
solution  in  powers  of  e must  fail  in  the  vicinity  of  the 
trailing  edge. 

* _ 

For  large  z it  is  also  shewn  in  Appendix  B that  F (/z,e) 
decays  like  (1/s  ) . Thus,  the  pressure  at  large  distances 

from  the  trailing  edge  is  given  by  the  lowest-order  result  with 
the  multiplicative  correi_tion  discussed  above  (see  (^.52)). 

The  Paradox  of  Infinite  Lift  - A Suggested  Experiment 

We  have  seen  that  in  the  absence  of  a shear  layer  the  total 
lift  on  any  section  of  the  trailing  edge  of  length  L grows  as 
v'L  . This  result  is  in  complete  agreement  with  the  exact  solution 
for  a finite  length  flapped  airfoil  (see  Ref.  1^,  e.g.).  One  can 
argue  that  in  reality  this  result  cannot  be  true.  The  total  lift 
must  either  reach  a maximum  or  decay  as  the  airfoil  chord  is  in- 
creased without  bound  for  fixed  flap  geometry.  The  flap 
eventually  becomes  totally  immersed  in  the  shear  layer  and  must 
Lecome  totally  ineffective  in  producing  lift.  Thus,  the  "steady 
ila'^ped  airfoil"  offers  a means  for  investigating  the  role  of  a 
i^hear  layer  in  limit i.ng  the  lift. 

The  paradox  could  perhaps  be  resolved  with  a very  simple 
experiment.  A two-dimensional  flat  plate  airfoil  could  be  fitted 
-..'ith  a flap  of  fixed  chord.  The  chord  of  the  airfoil  section 
forv;ard  of  the  flap  could  be  varied  while  the  total  steady  state 
lift  is  measured.  A plot  of  In(Lift)  versus  ln(c/£)  would  be 
expected  to  appear  something  like  that  shown  in  Fig.  9-  The  ratio 
c/£  \mere  the  lift  starts  to  depart  from  the  no-shear  result 
"straight  line  with  slope  1/2  in  Fig.  9)  would  be  a measure  of  the 
critical  rati-  of  shear  layer  thickness  to  flap  chord.  With  more 
ailed  theory,  this  point  could  also  be  calculated.  The  two-step 
r iel  does  not  appear  to  have  the  potential  for  describing  the 

..ration  because  the  velocity  at  the  surface  of  the  airfoil  can- 
* be  reduced  below  a fixed  value  as  the  shear  layer  thickness  is 
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increased.  A linear  velocity  profile  or  a power  law  profile 
would  be  more  realistic,  although  the  Wiener-Hopf  analysis  may 
become  intractable. 

The  Unsteady  Trailing  Edge 

To  continue  the  investigation  of  the  shear  layer  effect  on  the 
unsteady  trailing  edge  lift,  it  would  be  desirable  to  carry  out  the 
first-order  analysis  for  the  two-layer  model.  The  essential  and 
nontrivial  step  in  the  analysis  is  to  calculate  the  inverse  Fourier 

transform  of  the  shear  layer  function  Z(a)  given  by  (3.8)  and 

« 

(3»9).  Given  I (x)  the  first-order  pressure  is  a matter  of 
integration  (see  (3.35)  through  (3.38)).  These  results  would  be 
of  immense  value  in  assessing  the  role  of  the  shear  layer  on  the 
unsteady  lift. 


37 


V.  CONCLUSIONS 

We  have  Investigated  the  problem  of  an  oscillating  trailing 
edge  that  is  immersed  in  a two-dimensional  subsonic  shear  flow. 

The  ideal  problem  wherein  the  shear  layer  is  replaced  by  a two-step 
model  is  given  detailed  consideration.  An  approximate  solution  is 
developed  to  first  order  in  the  shear  layer  velocity  defect  (ratio 
of  displacement  thickness  to  boundary  layer  thickness)  and  is  valid 
for  arbitrary  length  scale  of  the  trailing  edge  deflection. 

The  general  solution  is  applied  to  calculate  the  steady  state 
lift  on  a linearly  deflected  flap  and  the  unsteady  lift  on  an 
exponential  flap.  The  latter  example  is  developed  to  first  order 
for  the  steady  state.  The  rpnciple  conclusions  of  our  analysis 
are  as  follows: 

1.  The  lowest-order  steady  state  lift  distribution  tends  to 
zero  as  the  square  root  of  the  distance  from  the  trailing 
edge  for  small  distance  and  as  the  reciprocal  square  root 
01  the  distance  for  large  distance.  The  linear  flap  has 
a weak  logarithmic  singularity  in  lift  at  the  hinge  line 
in  accordance  with  known  results. 

2.  The  total  lift  or.  a trailing  edge  section  of  length  L 
becomes  infinite  like  /L  as  L is  Increaved  without 

b bound.  This  result  is  in  agreement  with  tie  thin  six'foil 
solution  for  a finite  length  flapped  air full. 

3.  The  first-orcer  lift  on  the  steady  exponential  flap  is 
developed  in  detail.  It  is  shown  that  the  trailing  edge 
pressure  correction  varies  as  /jxj  and  is  attenuated 

by  a term  '.f  ffe  In  1/c)  where  e = 2S6/1  is  the  ratic 
of  shear  Isyer  thickness  to  flap  length.  We  ccnciude 
that  any  attempt  to  expand  the  solution  of  the  trailing 
edge  problem  in  a regular  perturbat icn  series  .mus*  fill 
near  the  trailing  edge. 
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Multiply  (A.l)  by  e and  integrate  over  all  x to  get 


«0 

1 f 

N,(x)  - ^ I 

— • X 


Since  (a  - is  analytic  in  the  lower  half  plane,  the 

function  N^(x)  is  zero  for  negative  x . For  x > 0 we  use 
the  contour  shown  in  the  above  sketch  to  evaluate  the  integral 


We  have 


Then 


\ix) 


r ‘ 

_ 1 I e^°*da 
^4  (a  - 


“1  * 


.2^iit/2 

t e 


ia.  x+iit/i»  r 

„ (,)  .is_i e-' 

1 ^ 


ia^x+in/4 


X > 0 


The  evaluation  of  NjCx)  is  similar  to  N^(x)  except  that 
the  contour  in  the  complex  o plane  must  be  closed  around  the 
branch  cut  from  into  the  lower-half  plane.  We  get 


NgCx) 


r 

i Trrrp^ 


-la,x-iT!/4 


X > 0 


, X < 0 


APPENDIX  B 


ASYMPTOTIC  ANALYSIS  OF 


F*(z,e) 


Consider  the  function 

P*(/z7e) 


with 


G(n,e) 


z 


e~^G(n»e) 


dn 


(B.l) 


(B.2) 


In  general,  P must  be  evaluated  numerically.  We  consider  below 
the  asymptotic  behavior  as  e -*  0 with  z « Oil)  . First,  we 
evaluate  G(n>e)  in  terms  of  known  functions.  We  have 


where 


G(n.e) 

Gq(ti,c) 

Gj^(n,e) 


- G^(n,e) 


(B.3) 


(B.il) 


(B.5) 


The  first  integral  is  easily  evaluated  in  terms  of  the  arc  tangent. 
We  get 

G^(n,e)  « 1 - tan”^  ~ (B.6) 

To  evaluate  the  second  integral  we  expand  the  integrand  by  partial 
fractions  and  translate  the  dummy  variable  of  Integration  in  each 
part . We  get 
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Gj^(n,e) 


«D 

-n 


-[e"^'Ej^(-il  - 


ie)  - 


00 

/ 


^(-n  + ie^ 


G(n,e)  “ 1 - I tan“^  | - G^(n,e) 

where  the  exponential  integral  is  defined  by 


larg  zi  < X , 


(B.7) 


(B.8) 


(B.9) 


To  obtain  the  asyniptotic  expansion  of  F for  ssrsali  e we 
must  first  develop  a uniformly  valid  expansion  of  G(n,e)  . We 
use  the  method  of  matched  asymptotic  expansion  to  accomplish  this 
task.  First,  we  expand  G(n»e)  in  the  limit  e 0 and 
*1  » OO)  . Then  we  expand  G for  c -*•  0 and  n ■ 0(e)  . The 
two  expansions  are  then  shown  to  have  a common  domain  of  validity 
by  matching.  Finally,  we  form  an  additive  composite  expansion 
with  the  two  results. 


Expand  G(nie)  for  e -»  0 


»o)g) 


We  use  known  properties  of  the  exponential  integral  to  cai‘ry 
out  the  following  program  (see  Ref,  15,  p.  226),  First,  we 


expand  S,  as  fellows: 
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Ej^(-n  i ie) 


•/ 


^ dt  + E^(-n  ± 10) 
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isi  r.£^ 

n J 1 - n 


Tr7?f  * ^1^-"  " 


to 


» -E,(n)  + iir  t le  I-  + O(e^) 
1 n 


(B.IO) 


The  last  result  is  valid  for  all  n greater  than  or  equal  0(1)  . 
Now  substitute  the  last  result  into  (B.7)  and  combine  the  result 
with  (B.3)  to  get 


Lin  G(n,e) 
£-*■0 


e 

ir 
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El(n) 
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+ O(e^) 


(B.ll) 


n>0(l) 


Expand  G(n«e)  for  e •*  0 , n = 0(e) 

For  n and  e both  small  we  can  write 

- ie)  *=  -Y  - Inl-n  ^ ie)  -n  s le  + O(e^)  (B.121 

where  y » .57721  is  Euler's  constant,  and 

ln(-n  1 ie)  = in(n“  + ± iv.  t ± tan~^  - (B.13) 

n 

We  use  (3.12)  and  (5.13)  in  (B.7)  and  (B.8)  to  get 

Liiu  G(n,e'  * ~ ^ 
e - 0 
n * Ole) 

which  result  is  regular  near  n = 0 but  is  of  0(e  In  c)  . 

Uniform  Expansion  of  G(n  ,e  ) 

Now  we  expand  (3.11)  for  small  n and  (B.14)  for  large  n . 

We  get 
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(B.14) 
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Llm  / Lxm  G ( n j e )\  = - — ( In  ti  + y - 1 ) 


(B.15) 


\ n>0(l)  / 

Lim  ,/Llm  G(n,£)\=  - ^(Inn  + Y-D 
ri-*0(l)'l  e-»-0  I 

\ n=  0(e)  / 


(B-.16) 


The  last  t.wo  results  are  identical  so  that  the  t.wo  asymptotic 
expansions  ;have  a common  domain  of  validity.  We  form  a uniformly 
valid  composite  expansion  by  adding  (B,ll)  and  (B. 14)  and  sub- 
tracting the  part  they  have  in  common  (either  (B.15)  or  (3.16)). 


The  final  result  is 


Lim  G( 
.e-»-0 

All  r. 


/ 2, 

it  ) 


(B.17) 


finally,  we  suost ituie  (B.17)  into  (B.l)  to  obtain 


Lim  r ( v/I,E)  =•■  + 


The  last  result  it  valuri  for  all  z and  in  general  must  be 
evaluated  r.umeracally  ^ Ur.iort unately , the  integral  is  not  sub- 
stantially more  Eimc.;c  than  the  original  integral  (B.l)  that  is 
valid  for  arbitrary  £ 

* 

V/e  conclude  cut  analysis  of  P with  results  for  small  and 
large  x ■ for  t - 0 , we  use  the  3 imit  value  of  tne  expression 
in  curly  braces  in  '’P  3 2).  We  get 
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The  functional  Tariation  ' s 

is  positive  and  the  amplitude 
the  contribution  to  the  integral  ; 
of  0(e/2^'^'^)  . Thus,  we  can  use 
in  the  approximate  form 


z - t The  function 
in  L . r't'r  a 3 ares  , 
, c <-t  near  n = z is 

5 . i b ' cr  t he  or 1 glnal  int  e gr a z 
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We  integrate  by  parts  in  (E.2C*  t 


P*(^.E)  - ^ 


Thus,  we  see  that 


i t'  z 

ft 


/-e*-:-.t,£-)  + KrlT?)  ' '(/s'?) 


an  i/*^  for  large  z in 

with  (B.irl 


decays 

contrast  to  Dawson’s  integral , We  sannct  evaluate  P 
since  we  have  already  neglected  terms  of  the  same  order. 





